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1  Introduction 


Ow  objective  is  to  describe  a  method  for  constructing  an  observer  for  the 
dynamical  system 

Ht)  =  /(^(0)  +  £ff*(*(0)“.(0.  x(0)=Xo,  (1) 

y(0  =  ^(2:(0). 

as  the  asymptotic  limit  of  nonlinear  filters  associated  with  the  “noisy”  ver¬ 
sion  of  (1)  :  '  '  I  ^  j  I 

di‘(t)  =  f{x*{t))dt  +  Yi9iix'it))ui{t)  +  ^/eN{x‘{t))dw{t),  (2) 

t=i 

a:‘(0)  =  xl, 

dC{t)  =  h{x^{t))dt  +  y/lRdvit),  e‘(0)=0, 

with  €  — >  0.  Here  x{t)  €  -K",  y{t)  €  Si’’  as  usual.  The  method  is  motivated 
by  some  large  deviation  results  of  Hijab  [3],  [4]  for  the  conditional  measures 
/-',5  of  (2). 

In  the  present  paper  we  present  results  of  this  general  method  as  applied 
to  the  linear  case  and  a  certain  class  of  nonlinear  systems.  The  general 
nonlinear  problem  will  be  treated  elsewhere. 

2  Observers  for  Linear  Systems 

In  this  section  we  provide  a  complete  description  of  the  method  as  it  applies 
to  linear  systems.  The  results  are  improvements  and  completions  of  earlier 
preliminary  accounts  provided  in  [l],  [2]. 

The  method  constructs  explicitly  an  observer  for  the  linear  system 

x[t)  =  Ax{t) x(0)  =  xo,  (3) 

v(t)  =  Cx{t), 

as  the  asymptotic  limit  of  (Kalman)  filters  for  a  family  of  associated  filtering 
problems 

di‘(0  =  Ax*[t]dt^-Bu[t)dt-\-yfiNdw{t),  x‘(0)=x*,  (4) 

de‘(0  =  Cx*{t)dt  +  y/lRdv{t),  C(0)  =  0. 
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Such  a  construction  is  suggested  by  the  fact  that  for  certain  choices  of 
Qq  =  cov(iq),  the  filters  are  independent  of  c,  as  discussed  in  Baras  and 
Krishnaprasad  [l].  Also,  the  solutions  of  (4)  converge  in  probability  as 
c  — »  0  to  the  solution  of  (3). 

The  work  of  Hijab  [3],  [4j  is  indispensible  here  in  deriving  a  large  de¬ 
viation  principle  for  the  conditional  measures  (see  Section  2.3),  and 
identifying  the  limit  of  the  filters  for  (4)  as  an  associated  deterministic 
estimator. 

2.1  Observers  and  Filters 

We  assume  as  usual  that  x{t)  €  .R",  u(t)  6  R"*,  y{t)  €  R^,  and  i  u(t) 
is  piecewise  continuous. 

Recall  that  the  observer  problem  consists  of  constructing  a  dynamical 
system 

rh{t)  =  Em{t)  4-  Fu[t)  +  Gy{t),  m(0)  =  mo,  (5) 

so  that  the  error 

e(t)  =  x(t)  -  Hm(t)  (6) 

decays  exponentially  fast  to  zero,  at  a  rate  controlled  by  the  designer,  in¬ 
dependent  from  the  choice  of  mo  and  Xo-  Here  the  matrices  E,F,G  and 
H  are  possibly  time-varying  and  the  dimension  of  m(t)  is  not  necessarily 
n.  This  of  course  reflects  the  fact  that  the  initial  condition  Xo  is  unknown, 
and  the  best  that  can  be  done  is  to  approximately  estimate  i(t)  by  Hm[t) 
in  this  way. 

Solutions  to  this  problem  are  well  known,  first  given  by  Luenberger  [6]. 
In  particular,  if  the  pair  {C,A)  is  detectable,  then  there  exists  a  matrix  T 
such  that  the  matrix  A +  TC  has  eigenvalues  in  the  open  left  half  plane. 
Then  set 

E  =  A  +  TC,  F  =  B,  G=-T,  H  =  I. 

In  this  case  the  error  (6)  satisfies 

«(t)  =  (A -f  rC)e(t),  e{0)  =  Xo  -  2o, 

and  the  eigenvalues  of  A  4-  TC  can  be  arbitrarily  assigned  by  the  designer 
if  and  only  if  (C,A)  is  observable. 
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Consider  the  system  (3).  Define  f(t)  =  /o*y(«)ds,  so  that  (3)  becomes 

i(t)  =  Ax{t) ->r  Bu[t),  x(0)  =  Xo,  (7) 

i{i)  =  Cx(t),  e(0)=0. 

Then  associate  with  (7)  the  family  of  filtering  problems  (4),  where  w,  v  are 
independent  standard  A:-dimensional,  respectively  p-dimensional  Brownian 
motions.  The  initial  condition  xj,  is  Gaussian,  independent  from  v,  w  with 
E{Xq)  =  /Iq,  cov(xo)=Qo  ,  where  Qq  is  positive  definite.  Note  that  the 
(small)  parameter  c  controls  the  intensity  of  the  noise.  The  matrix  R  is 
assumed  positive  definite. 

As  is  well  known,  the  minimum  variance  estimate  i‘(t)  =  E(ilt)  | 
^  ^  0  linear  Gaussian  filtering  problem  (4)  is  given  by 

the  Kalman  filter  [7] 

dx‘(0  =  Ar(t}dt  +  Bu(t)dt  +  Q‘(t)C'(RR’)-^C(d^*(t)  -  Cx^(t)dt), 

i-(0)  =  (8) 

where  Q'  satisfies  the  Riccati  equation 

(3‘(0  =  AQ^(t)^Q‘(t)A’-Q‘(t)C'(RR')-^CQ^(t)  +  NN\  (9) 

Q‘(0)  =  Ql/e. 

Note  that  these  filters  depend  on  c  only  via  the  matrix  Ql/e.  In  fact,  if 
we  choose  Qq  =  eQo,  then  ail  the  filters  are  independent  of  c  and  identical 
with  the  filter  for  f  =  1. 

Following  Hijab  [3],  it  is  convenient  to  consider  the  filter  (8),  (9)  as  a 
map 

e{«),0<a<t  8-^  x‘(0. 

2.2  Deterministic  Estimation 

Following  Mortensen  [8]  and  Hijab  [3],  we  associate  with  (7)  the  determin¬ 
istic  (noisy)  system 

i{t)  =  Az{t)  +  Bu(t)  +  Nw{t),  s(0)  =  xoi 
f(0  =  Cz{t)  +  Rv{t),  f(0)=0. 


(10) 


and  energy  cost  functional 


1  1 

J,(zo,w,v)  =  -(^0-  u)  Qo  (*0  -  m)  +  2  y  (“'(«)'*"(«)  +  <'(«)'v(5))  ds, 

.  /  . 

where  t  •— »  w(t)  £  JR*,  t  •-*  t>(f)  €  are  piecewise  continuous,  the  rank  of 
N  is  n,  and  Qo  is  positive  definite. 

A  minimum  energy  input  triple  (ao,ti;*,t;*)  given  f(s),  0  <  5  <  t,  is 
a  triple  that  minimises  J(  subject  to  (10)  and  produces  the  given  output 
record  f(s),  0  <  s  <  <.  The  deterministic  or  minimum  energy  estimate  of 
z(t)  given  f(s),  0  <  s  <  t,  is  the  endpoint  z(t)  of  the  trajectory  2*(5),  0  < 
s  <  t,  of  (10)  corresponding  to  a  minimum  energy  input  triple:  z(t)  =  z*(t). 
According  to  Krener  [9],  i  is  the  solution  of  the  Kalman  filter  equations 

i(0  =  Ai(0  +  Bu(0  +  g(0C'(RRrMf(0-Ci(0),  (12) 

i(0)  =  M. 

Q(t)  =  AQ(t)  +  Q(t)A'-Q(t)C'(JiJi')-^CQ(t)-hNN',  (13) 

Q(0)  =  Qo. 

As  in  the  stochastic  case  (Section  2.1),  it  is  convenient  to  consider  the 
deterministic  filter  (12),  (13)  as  a  map 

/:C>(|0,t],«'’)  — 
f(s),  0  <  a  <  f  I — ►  i(t). 

Note  that  the  deterministic  filter  coincides  with  the  stochastic  filter  for 
<  =  1,  that  is,  7*.  Also,  z(t)  is  obtained  from  an  optimal  control  problem, 
and  is  determined  by  a  Hamilton-Jacobi-Bellman  equation  [3],  [8]. 

We  now  prove  that  as  e  — *  0  the  stochastic  filter  7'  (8),  (9)  converges 
to  the  deterministic  filter  7  (12),  (13). 


Theorem  1  Suppose  that  (4)  has  initial  conditions  zj,  Gaussian  with  mean 
/i'  and  covariance  Q\,  satisfying 


nrn-<?‘  =<?o. 

*-•0  ( 


where  Qo  is  positive  definite.  Then 

iimf:|z*(0-i(0  |*=o 

uniformly  in  t  £  [0,7’]. 


Proof:  Now  Q‘{t)  — »  Q(/)  uniformly  in  t  G  [0,  Tj.  Applying  Ito’s  rule  to 
I  f'(<)  —  2{t)  I*  and  taking  expectations,  we  find  that,  given  6  >  0, 

E  I  x^it)  -  z{t)  I’  <  (I  /x*  -  /X  |»  +6K  (1+  I  /x*  1)  +  eK)  e^\ 

for  all  e  sufficiently  small,  where  K  >  0.  The  desired  result  follows  from 
this  inequality.  □ 

2.3  Large  Deviations 

Consider  the  stochastic  differential  equation  (4),  with  initial  condition  Zq  = 
zo  for  all  «  >  0.  In  this  section  we  take  u  =  0.  Let  P‘  be  the  probability 
measure  induced  on  fl"  =  C([0, Tj,!?")  by  the  diffusion  z*.  It  is  well  known 
from  the  theory  of  Ventcel-Friedlin  (see  Varadhan  [10])  that  the  family  of 
measures  satisfy  a  large  deviation  principle.  Moreover,  as  £  — »  0,  P* 
converges  weakly  to  the  degenerate  measure  concentrated  on  the  unique 
solution  z  of  (3). 

We  now  consider  the  observation  equation  in  (4).  Let 
unnormalised  conditional  measure  on  fl"  of  z*  given  ^  €  fl’’  =  C((0,r],]R'’) 
where  the  diffusions  are  initialised  as  above.  For  a  “control”  t  — ►  tx)(f), 
let  Zfu  denote  the  unique  solution  to  (10).  The  function  t;  is  defined  by 
-^”*(^(0  ~  C'e„;(t))  when  {  is  C*.  Hijab  [4]  proved  the  following  large 
deviation  result  for 

Theorem  2  For  any  open  subset  0  and  any  closed  subset  C  of  H", 
li^inf  «logQ;|((,,p)(0)  >  -I  (zo,  f,  0) 

limsupelog9;i((,„)(C)  <  -/(xo.^C) 


where  for  A  C  fl", 


nxo,i,A) 


(u;(s)'ty(5)  +  Z^{syC'Cz^{s))  ds  (14) 


=  inf  1 5/’ 

«"  (2/0 

-  z^{syc'd^{s)  I  «„,(0)  =  xo,  z„;  e  >l| 


Proof:  Define,  for  each  {  €  n’’,a;  G  fl", 

4>iu,0  =  -HTyCijiT)  + 

(^CityCAu;{t)  +  K;{tyC'Cu;[t)  -  \i{ty C N N' C i{t))  dt. 

There  exist  constants  A,B  depending  only  on  such  that 

<  >4  +  -B||w||. 

Then  arguing  as  in  Varadhan  [11], 

lim  lim  sup  flog  f  exp  f— =  ~®o. 

R^oo  /{«. -4(«,f)>/e}  \  t  ^  'V  * 

But  this  estimate  is  enough  to  prove  the  theorem.  See  Hijab  [4]  and  Varad- 
han  [11]  for  details.  □ 

The  minimisation  in  (14)  is  an  optimal  control  problem  similar  to  the 
one  in  Section  2.2,  but  with  fixed  initial  condition  zq.  James  and  Baras  [13} 
have  made  a  simple  generalisation  of  Theorem  2  in  which  the  variational 
problem  arising  is  exactly  the  optimal  control  problem  for  deterministic 
estimation  in  Section  2.2. 

Assume  that  the  initial  conditions  Zq  of  (4)  have  (unnormalised)  density 
5‘(zo)  =  exp^-^(zo-M)'Co*(^o-Ai)). 

Let  unnormalised  joint  conditional  measure  of  (3:*,Zo)  on 

fl"  X  J?"  given  f  €  fl'’.  The  following  result  is  quoted  from  (13). 


Theorem  3  For  any  open  subset  0  and  any  closed  subset  C  of  fl”,  and 
for  any  open  subset  Cq  and  any  closed  bounded  subset  Co  of  JR'^,  we  have 

lim  mf  e  log  {0  x  Oo)  >  -J  {0  x  Oo,  0 

£  log  (C  X  Co)  <  -J{C  X  Co,  0 

where  for  A  X  Ao  C  X  JR**, 

JiAxAo,C)  =  Jnf^|i(xo-/i)'Qo*(^o-M)  +  (15) 


This  theorem  implies  that  if  ^  f  is  an  actual  output  record  of  the 
system  (10),  then  as  e  — ♦  0,  converges  weakly  to  a  degenerate  mea¬ 

sure  concentrated  on  the  optimal  initial  condition  and  optimal  trajectory 
0  <  s  <  T  of  (10)  corresponding  to  a  minimum  energy  input  triple. 
As  pointed  out  in  Section  2.2,  the  deterministic  estimate  of  the  state  at  time 
T  is  a  functional  of  this  optimal  path,  namely  its  value  at  T:  z{T)  =  z*{T). 

Thus  in  a  weak  sense,  x*(T)  -+  z{T),  and  the  large  deviation  principle 
fo*’  ^(*.*o)|f  characterises  the  limiting  filter  as  the  deterministic  filter. 

2.4  Observer  Design 

From  Sections  2.2  and  2.3  it  is  plain  that  the  deterministic  estimator  (12), 
(13)  is  a  natural  candidate  for  an  observer  for  the  linear  system  (3).  We 
make  the  natural  assumption  that  the  pair  (C,  A)  is  detectable.  Recall  that 
N  has  rank  n  and  R  is  positive  definite.  The  design  parameters  are  Qo,  N,  R 
and  p,. 

Then  from  (12)  we  define 

m(t)  =  Am(t)-fBu(t)4-Q(t)C'(RR')-‘(y(t)-Cm(t)),  (16) 

m(0)  =  mo  =  p, 

where  Q(t)  is  the  solution  of  the  Riccati  equation  (13).  The  inverse  P(t)  of 
Q(f)  is  the  solution  of 

P(t)  =  -P(t)A-A'P(t)-P(t)NN'P(t)  +  C\RRy^C,  (17) 
m  =  P„  =  q;‘. 


Since  we  are  interested  in  the  asymptotic  behaviour  of  e(f)  =  x(<)  — m(l), 
it  is  important  to  obtain  bounds  for  ||  Q[t)  ||,  ||  P{t)  1|.  To  this  end  we 
interpret  Q(t),  P(t)  in  terms  of  control  problems.  Write  H  =  R~^C. 
Consider  the  control  problem 

-  ri  =  A*t)  +  H'v,  T){T)  =  h,  (18) 

where  h  is  given  and  v  is  the  control.  We  minimise 

Ji(v)  =  f?(0)'(?on{0)  +  f  (v{t)'v(f)  +  r){tyNN'T){t))  dt.  (19) 

Jo 

Then  the  optimal  control  for  (18),  (19)  is  given  by  the  following  algorithm. 
Consider  the  system  of  equations: 


A  =  AX  +  NN’fi,  A(0)  =  Qo^(0), 

-h  -  A'fi~  H’HX,  t){T)  =  h. 

Then  an  optimal  control  is  v(t)  =  Moreover, 

min  Ji(v)  =  h'Q(T)h  =  X(T)h. 

In  addition,  the  following  relation  holds: 

=  Q(0’7(0»  for  all  t, 

where  (?(t)  is  the  solution  of  the  Riccati  equation  (13). 

Similarly,  consider  the  control  problem 

A  =  ^A  +  7Vv,  A(r)  =  /i. 

Again  h  is  given  and  v  is  the  control.  We  minimise 

J,(t;)  =  A(0)'PoA(0)  +  r  (v(l)'v(l)  +  A(t)'H'HA(t))  dt. 

Jo 

The  system  of  necessary  conditions  is  given  by 

A  =  AA  +  NJV'r),  A(r)  =  h, 

-  ^  -  n'HK  17(0)  =  PoA(O), 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


and  an  optimal  control  is  v(t)  =  N'r}(t),  with 

min  J2(t;)  =  h'P(T)h  =  rj[T)h,  (26) 

and 

^(t)  =  P(t)X(i),  for  all  i,  (27) 

where  P(t}  is  the  solution  of  (17). 

Since  P  is  positive  definite,  in  particular  nonsingular,  the  pair  (P,A)  is 
detectable.  Thus  there  exists  a  matrix  A  such  that 

rj'(A  +  AH)T}  <  -ao  I  »7  I*,  ao  >  0.  (28) 

Also,  since  A'  has  rank  n,  the  pair  (A,N)  is  controllable,  and  there  exists 
a  matrix  F  such  that 

A'(A  +  A^r)A  >  /?o|Al’,  /?o>0.  (29) 

Theorem  4  Under  the  above  assumptions,  rvt  have: 

II  QiT)  II  <  (ll  Qo  il  +  (30) 

II  P{T)  II  <  (||  Po  II  +  -  — )  =  P 

Proof:  Consider  in  (18)  a  feedback  control 

v{t)  =  A'r?(t). 

The  corresponding  state  is  the  solution  of 

-q  =  {A’  +  H'A')rj,  v{T)  =  h.  (32) 

Therefore 

h'Q{T)h  <  rj(0)'Qo»7(0)  +  T  »?(0'  ri[i)dt.  (33) 

JO 
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2.5  Convergence  of  the  Linear  Observer 

We  now  use  the  bounds  (30),  (31)  to  prove  the  following. 


Theorem  5  The  dynamical  system  (16),  (13)  is  an  observer  for  the 
linear  control  system  (3)  provided  that  [C,A)  is  detectable  and  the  above 
assumptions  hold.  That  is,  there  exists  constants  K  >  0,  7  >  0  such  that 

I  x{t)  -  m{t)  I  <  K  I  Xo  -  mo  I  e"'’* 


for  all  t  >  0. 

Proof:  From  (30),  (31)  it  follows  that 

1  P(‘)A  1  > 

(36) 

and  ^ 

<  yP(t)X  <  plAl^ 

9 

(37) 

Now  e(<)  =  x{t)  -  m(()  satisfies 

e{t)  =  iA-Q{t)H'H)e{t). 

Using  (17),  (36)  we  deduce 

=  -c[t)'(P{f)NN'P{t)+H'H)c{i) 
<  -e{tyP[t)NN'P[t)e{t) 


where  NN'  >  tqJ,  tq  >  0.  This  together  with  (37)  implies 

at  pq* 

Set  7  =  ro/2p5*.  Therefore 

e(tyP(t)e(t)  <  e(0)'Poe(0)e-*'»‘ 


S  ‘ii'i 


and 


I  e(0  |»  <  gc(0)'Poc(0)c-*"‘ 

from  which  the  desired  result  follows.  □ 

Finally,  we  state  the  following  result  which  is  a  consequece  of  standard 
facts  concerning  the  Riccati  equation  (13). 

Theorem  6  Given  the  linear  system  (3),  where  {C^A)  is  detectable,  an 
n  X  m  matrix  N  such  that  {A,N)  is  stabilisable,  and  a  positive  definite 
matrix  R,  then  there  exists  a  unique  non-negative  definite  solution  Q  to 
the  algebraic  Riccati  equation 

AQ  -i-  QA'  -  QC'{RR')-^CQ  -h  NN'  =  0, 

the  matrix  A  —  QC'{RR')~^C  is  exponentially  stable,  and  the  system 

m{t)  =  Am{t)-\-Bu{t)-^QC'{RR')-^iy{t)-Cm{t)), 

m(0)  =  mo, 

is  a  time-invariant  observer  for  the  given  system. 

3  Observers  for  Nonlinear  Systems 

We  consider  a  nonlinear  dynamical  system  with  linear  observations: 

i(0  =  /(^(O).  ^(0)  =  ^0,  (38) 

V(0  =  C’a:(0. 

We  assume  that  /  :  i?”  — »  M”  is  smooth  with  bounded  derivatives,  and 
write 

A{x)  =  Df[x) 
for  the  n  X  n  matrix  of  first  derivatives.  Set 


3.1  Observer  Design 

Motivated  by  the  linear  design,  we  construct  an  observer  for  (38)  as  an  ap¬ 
proximation  to  the  corresponding  deterministic  estimator.  Associate  with 
(38)  the  system 


^(0  =  /(^(O)  +  •/Vu;(t),  a(0)  =  Zo,  (39) 

C(t)  =  Cz(t)  +  Jiv(t),  f(0)=0, 

where  rankA^  =  n,  H  is  positive  definite,  and  the  energy  functional 
1  If* 

Jt{zQ,w,v)  =  -{zo-tiyPoizQ~(i)  +  -  /  (u;(s)'u;(s) -t- v(s)'v(5))  ds,,  (40) 

Z  Z  JO 

where  Po  is  positive  definite. 

According  to  Hijab  [3],  the  deterministic  estimate  z  is  the  solution  of 

i(<)  =  /(f(<))  +  e(‘)c'(Mr‘ (r(0  -  C»(1)) ,  (41) 

i(0)  =  <J. 


where 

Qit)-^  =  D^S(z(t),t), 

and  5(2,  t)  is  the  solution  of  the  Hamilton- Jacobi-  Bellman  equation 

^^Siz,t)  +  ff(z,t,DS(z,t))  =  0,  (42) 

5(2,0)  =  ^(«o  -  m)'Po(*o -m), 

where 

ff(z,t,a)  =  a/(z)  -I-  ^aNN’a'  -  lz'C'{RR')-'Cz  -f  f(t)'C2. 

In  the  linear  case  the  solution  of  (42)  is  a  quadratic  form  and  Q{t)  = 
P(t)~^  satisfies  a  Riccati  equation.  However,  in  the  general  nonlinear  case, 
solutions  are  not  smooth  and  must  be  interpreted  in  the  viscosity  sense. 
Thus  (41)  is  not  well  defined  in  the  large.  We  seek  therefore  to  approximate 
S(z,t)  by  a  quadratic  form,  and  replace  the  Hamilton- Jacobi  equation  (42) 
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by  a  simpler  Riccati  equation.  In  this  way  we  will  obtain  a  well  defined 
observer.  Write  H  =  R'^C. 

Suppose  that  S  is  smooth  in  a  neighbourhood  of  (MiO).  Denoting  com¬ 
ponents  by  superscripts  and  partial  derivatives  by  subscripts,  and  using  the 
summation  convention,  for  small  t  we  have  at  (l(t),t): 

-S,i{NNYSi,  + 

using  the  fact  that  5t(2(t),t)  =  0  from  the  definition  of  i(t).  If  5  were 
quadratic,  the  third  order  terms  vanish.  This  suggests  that  we  replace  (41) 
by 

mit)  =  /(m(t))-f-Q(t)C'(RR’rU!/(t)-Om(t)),  (43) 

m(0)  =  mo  =  M. 

where  now  Q(t)  =  and  P(t)  satisfies  the  Riccati  equation 

P(t)  =  -P(t)A(m(t))  -  A(m(t)yP(t)  ~  P(t)NN'P(t) -h  ff'R,  (44) 
P(0)  =  Po. 

Also  Q(t)  is  the  solution  of 

Q(t)  =  A(m(t})'Q{t)-i-Q(t)A(m(t))-Qit)H'HQit)  +  NN\  (45) 

g(0)  =  Qo  =  Po'- 

Once  again  it  is  important  to  obtain  bounds  for  ||  Q{t)  ||,  ||  P(t)  ||.  To 
recover  estimates  similar  to  (30),  (31),  we  assume  that  the  pair  (II,A(x)) 
is  uniformly  detectable,  that  is,  there  exists  a  bounded  Borel  matrix  valued 
function  A(i)  such  that 

n'(A(i) -i- A(i)H)t}  <  -QoItjI*,  Qo  >  0,  (46) 

for  all  X  €  Ft".  Since  N  has  rank  n  and  ||  A  ||  <  oo,  the  pair  (A(i),A') 
is  uniformly  controllable,  and  thus  there  exists  a  bounded  Borel  r(i)  such 
that 

X'  {A(x)  +  Arr(x))  A  >  ^0  I  A  I*,  00  >  0,  (47) 
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■V 


j',-.  it".V .  A>i 


sup 

«€#t" 


r(x)  II . 


for  all  I  6  JR''.  Set 


II  A  II 


sup  II  A(i)  II, 


r  II  = 


Then  using  the  methods  of  Section  2.4,  the  following  generalisation  of 
Theorem  4  can  be  proven. 


Theorem  7  Under  the  above  assumptions,  we  have: 


II  Q{T)  II  < 

(ll  Oo  II 

,  l|N|l*  +  ll  Al|»\ 

20.  j  -  ’■ 

(48) 

li  P{T)  II  < 

(ii  n  II 

,  II  »  P  +  II  r  ll’'! 

^  2^0  J  - 

(49) 

3.2  Asymptotic  Convergence 

We  wish  to  prove  that  the  system  (43),  (45)  is  an  observer  for  the  nonlinear 
system  (38).  This  is  possible  provided  that  we  can  bound  the  region  where 
the  initial  condition  lies  and  provided  the  second  derivative  of  /  is  not  too 
large. 

Consider  £>v4(i)  =  r>*/(z).  For  any  x,  I?*/(i)  £  L  (JR", L  (iR",  2R")) 
and  we  denote  ||  D^f  ||  the  suprimum  over  i  of  the  norm  of  the  linear 
operator  D^f{x). 

Note  that  the  numbers  p  and  q  defined  by  (48),  (49)  are  functions  of 
the  design  parameters  Pq,  N,  and  R;  and  the  given  data  /  and  C.  The 
designer  is  free  to  chose  the  parameters  within  the  stated  constraints.  Also 
NN'  >  tq/  for  some  Tq  >  0.  Define 


r-m.t 


<p{Po,N,R)  = 


To 


rp 


pUi 

*r\ 


Theorem  8  Assume  that 


To  -  mo  III  D*/  II  <  <p{Po,N,R). 
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(50) 


(51) 


Then  the  dynamical  system  (43),  (45)  is  an  observer  for  the  nonlinear 
system  (38)  provided  that  (C, >4(i))  is  uniformly  detectable  and  the  above 
assumptions  hold.  That  is,  there  exists  constants  K  >  0,  “7  >  0  such  that 


x(t]  -  m(t)  I  <  JC  I  xd  -  mo  I  e 


for  all  t  >  0. 


Proof:  Now  e{t)  =  x{t)  —  m(t)  satisfies 

e(t)  =  f(x(t))-f(m(t))-Q(t)ff'ffe(t). 

From  (44)  we  deduce 

^,(typ(t)e(t)  =  -e(0'  (P(t)NN>P(t)  +  H'ff)  e(t)  (52) 

+2e(0'P(<)  ifix(t))  -  f(m(t))  -  Q(t)H'ffe(t)) 

=  -e(ty  (P(t)NN'P(t)  +  ff'ff)  e(t) 

■+2e{tyP{t)  f  f  rD^f  {m{t)  +  rsx{t))  e{tydrds 

Jq  Jq 

<  ,(()'  (-p  +  I  ^”'’(0^(1)  I  li  CV  ll)  <(l).  (53) 

By  the  assumption  (51)  we  can  find  Po,N,R  such  that 

I  e(0)  III  i?V  II  <  <p[Po,N,R), 

hence 

1  p„‘'Mo)  III  DV II  < 

or 

+  I  I  II -dV  II  <  0. 
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Since  P’/*(()e(t)  is  continuous,  there  exists  an  interval  [0,<o)  such  that 

-  T  +  I  I  p'"  II  £>■/  II  <  0, 
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on  [0,  to)-  But  from  (53), 


±  I  P^'Ht)e(t)  I  <  0 


on  (0,<o)i  and  thus 


I  I  <  I  Po'''^(o)  I 

on  fO,  to)  .  By  continuity  we  have 

I  P''’(lo)«(<o)  I  <  I  Fo’'’<(0)  I, 

and  we  can  proceed  from  to  on.  Therefore  in  fact 

I  I  <  (p  - ‘ 

for  all  t  >  0,  and  (53)  implies 

<  -{ I  «(i)  1= . 

But  from  (49), 

c(t)'P(t)e(t)  <  II  P(0  III  c(0  I*  <  P|e(0r. 

hence 

which  implies 

e(t)'P(t)e(t)  <  e(0)'Poc(0)c"^. 

Therefore,  using  (48), 

i«(oi*  <  II  c?(o  II 

<  ge(t)'P(t)e(t) 

<  ge(0)'Poe(0)e~^‘, 

from  which  we  deduce  the  desired  result.  □ 

Remark  It  follows  from  (52)  that  if  the  design  parameters  Po,N,R  were 
chosen  so  that 

0  <  «  = 

then  the  assumption  (51)  is  unnecessary.  Then  (43),  (45)  is  an  observer  for 
(38)  independent  of  the  initial  conditions.  Unfortunately  this  inequality  is 
at  best  difficult  to  achieve. 
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